Elliptic Quantum Billiard 



Holger Waalkens, Jan Wiersig, Holger R. Dullin 

Institut fur Theoretische Physik and 
Institut fur Dynamische Systeme 
University of Bremen 
Postfach 330 440 
28334 Bremen, Germany 

E-mail: jwiersig@physik.uni-bremen.de 

February 5, 2008 



Abstract 

The exact and semiclassical quantum mechanics of the elliptic billiard is investigated. The classical 
system is integrable and exhibits a separatrix, dividing the phase space into regions of oscillatory 
and rotational motion. The classical separability carries over to quantum mechanics, and the 
Schrodinger equation is shown to be equivalent to the spheroidal wave equation. The quantum 
eigenvalues show a clear pattern when transformed into the classical action space. The implication 
of the separatrix on the wave functions is illustrated. A uniform WKB quantization taking into 
account complex orbits is shown to be adequate for the semiclassical quantization in the presence of 
a separatrix. The pattern of states in classical action space is nicely explained by this quantization 
procedure. We extract an effective Maslov phase varying smoothly on the energy surface, which is 
used to modify the Berry- Tabor trace formula, resulting in a summation over non-periodic orbits. 
This modified trace formula produces the correct number of states, even close to the separatrix. 
The Fourier transform of the density of states is explained in terms of classical orbits, and the 
amplitude and form of the different kinds of peaks is analytically calculated. 



1 Introduction 



The semiclassical quantization of a Hamiltonian system is deeply connected to the structure of its 
phase space. The generic Hamiltonian system contains a complicated mixture of near-integrable and 
chaotic motion, and a consistent semiclassical quantization scheme does not exist for the generic 
case. It does exist, however, in the non-generic limiting cases, the integrable and the ergodic 
systems. 

During the last two decades much progress has been made for ergodic systems. The Gutzwiller 
trace formula gives the density of states as a sum over periodic orbits [jl], Q. It works well if all 
periodic orbits are unstable and isolated, i.e. for hyperbolic systems, and is the starting point for 
most calculations in this field. 

The investigation of integrable systems reaches back to the beginning of quantum mechanics. 
Bohr and Sommerfeld, among others, succeeded in the quantization of actions. The hydrogen 
atom is the most famous example. The deficiency of the old quantum mechanics to calculate, 
e.g., the spectrum of the helium atom was set into clear light by Einstein Q. He formulated the 
Bohr-Sommerfeld quantization conditions in terms of invariant tori which foliate the phase space 
of integrable systems. Moreover he noted that this foliation is absent in generic systems, such that 
this quantization scheme fails. The same year when Schrodinger introduced his famous equation 
Brillouin explained that the quantization of tori is a consequence of the single-valuedness of the 
quantum mechanical wave functions. Keller showed how the quantization conditions have to be 
modified because of caustics in the classical motion 0. For an integrable system with / degrees of 
freedom described by its phase space variables (q, p) = (gi, qf,pi, ■■■,Pf) the Einstein-Brillouin- 
Keller (EBK) quantization conditions for the actions Ij read 

If ai 

where the integration is to be taken along / topologically independent paths ji around the /-torus, 
and ai are the Maslov indices due to the classical caustics. 

This paper deals with a nontrivial example of an integrable system: the planar elliptic billiard. 
Classically "billiard" refers to a system in which a point particle moves freely inside a domain and 
is elastically reflected at the domain boundary. The corresponding quantum mechanical problem is 
to solve the Laplace equation inside the domain with Dirichlet boundary condition. Billiards have 
become popular because for them on the one hand the classical calculations are easier than for 
systems with smooth potential, and on the other hand they allow for experimental measurements, 
e.g. H 0, ||. Jacobi || showed that elliptic coordinates separate the geodesic flow on the ellipsoid, 
which contains our billiard as a limiting case. The same coordinate system leads to the separation 
of Schrodinger's equation, including the Dirichlet boundary condition, into two Mathieu equations. 

Integrability allows for a semiclassical quantization a la EBK. In our case, however, there are 
problems due to the presence of a separatrix. For all energies it divides phase space into regions of 



rotational and oscillatory motion [ 10 ] . The Maslov indices are different for the two classical regions, 
resulting in a discontinuity in the EBK quantization condition which in turn leads to ambiguities 
for states close to the separatrix. To overcome this problem it is necessary to introduce a uniform 
quantization condition. This can be done by investigating the asymptotics of the solutions of 



Schrodinger's equation close to the separatrix, which was carried out for the elliptic billiard in [11|. 



We will follow a different approach based on connection matrices between the amplitudes of WKB 



wave functions, see [12, 13|, the review of Berry and Mount [14], and the references therein. This 
approach can be interpreted in terms of orbits with complex classical action. A transformation of 
the quantum mechanical eigenvalues to classical action space allows for a unique mapping of the 
eigenvalues to the quantum numbers. 

The discontinuity in the EBK quantization condition carries over to the Berry- Tabor trace 
formula [|D| - the analogue of the Gutzwiller trace formula for integrable systems. The resonant 
tori, foliated by families of periodic orbits, take over the role of the isolated periodic orbits as the 



objects to be summed over. The uniformization employed in [15] incorporates the influence of orbits 
with negative classical action in the vicinity of isolated stable periodic orbits. We will introduce 
a modification of the Berry- Tabor trace formula which takes care of the separatrix in terms of an 
effective Maslov phase varying smoothly across. It will turn out that in order to correctly produce 
eigenstates close to the separatrix the summation has to be taken over non-periodic orbits, i.e. over 
classical tori with in general irrational winding number. This trace formula incorporates all kinds 
of classically non-real orbits. 

The Berry- Tabor trace formula was further investigated by Richens fl6(| , who showed that it 
contains contributions of the stable isolated periodic orbits, whose contributions are equal to the 
corresponding terms in the Gutzwiller's trace formula §§. We will use his results in the study of 



the length spectrum, extending the "inverse quantum chaology", see e.g. [17, 18], to integrable 
systems. 



After completion of this work there independently appeared a preprint by M. Sieber [19], whose 



first part contains considerations similar to parts of our paper. Nevertheless, the key point of his 
paper being the semiclassical consequences of the deformation of an elliptic billiard to an oval, 
while we concentrate on the separatrix and on complex orbits, the overlap is only mild. 

The organization of our paper is as follows. We start with a short summary of the classical facts 
in sec. || and perform the exact quantum mechanical calculations in sec. ||. In sec. |] we introduce a 
uniform WKB quantization condition. This method gives an effective Maslov phase which is used 
to modify the Berry- Tabor trace formula for systems with a separatrix in sec. ||, resulting in a 
sum over non-periodic orbits. The appearance of resonant tori in the length spectrum, i.e. in the 
Fourier transform of the density of states, is investigated in sec. |6[ The conclusion and an outlook 
are given in sec. |?]. 

2 Classical mechanics 

The classical dynamics of the planar elliptic billiard has been investigated by many authors, e.g. [p0| 



10, 21, |22|] and the references therein. We just give a summary of the facts important for our 



purpose. Scaling the longer semimajor axis to one, the boundary of the billiard is described by 

v 2 

x 2 + = 1, < a < 1 , (2) 

1 — cr 

with foci at (x,y) = (±a,0). The boundary of the ellipse is a ^-coordinate line of the elliptic 
coordinates (p, (p) given by 

(x,y) = (a cos cosh />, asin^sinhp) , (3) 

where the coordinate ranges are 

< p < Pmax = arccosh(l/a), < (J) < 2ir . (4) 

The lines p = const are confocal ellipses and the lines <ft = const are confocal hyperbolas. Intro- 
ducing the conjugate momenta (pp,p<^), a reflection at the boundary p = Pmax is simply described 
by (p, — ► (p, tfi, —p p ,P(j,)- The Hamiltonian of a freely moving particle of unit mass reads 
in these coordinates 

H - P ? +P t (5) 

2a 2 (cosh 2 p — cos 2 cf>) 

Multiplying Eq. (|5l) with the denominator of the right hand side yields the separation constant K: 



p 2 + (K- 2Ea 2 cosh 2 p) =0, 
V$ ~{K - 2Ea 2 cos 2 (f>) = . 



(6) 
(7) 




Figure 1: Effective potentials V e s(p) and V e ff(</0 and the tori of the elliptic billiard. Real tori are marked by R 
(rotations) or O (oscillations). The classically forbidden motion are TR (tunneling from region 1 to 2 in the radial 
direction), TO (tunneling from region 1 to 2 in the angular direction), and O c and TO c , the complex continuations of 
O and TO below the minimum of V c s(<j>). Orbits with complex position are indicated by a circled cross. The reason 
for taking negative values of p into account will become clear in sec. El 



Here E > is the energy, and k 2 = Kj (2E) is the second constant of the motion. Since H and K 
are in involution, the system is integrable and therefore the energy surface is foliated by invariant 
Liouville 2-tori. Each of the equations (||) and (^) can be interpreted as a Hamiltonian system 
with one degree of freedom, with effective energy E p = —Ek 2 and Em = En 2 respectively, and a 
sum of kinetic term and effective potential V e g(p) = — Ea 2 cosh 2 (p) or V e s((f>) = Ea 2 cos 2 (4>). The 
topologically different types of motions can be discussed in terms of the effective potentials shown 
in Fig. [p. It is obvious that there are only two types of classically allowed motion. For k 2 > a 2 the 
trajectories avoid the interior of the ellipse cosh(p) = k/cl, touching its boundary between every two 
consecutive reflections at the billiard boundary p = /0 max - This "type R" motion is similar to the 
rotational motion in a planar circular billiard. For n 2 < a 2 the trajectories always cross the x-axis 
between the foci; they are confined to the domain enclosed by the hyperbolas cos (ft = ±k/ci. This 
"type O" motion involves an oscillation in (ft. The special values k? = and k 2 = 1 represent the 
stable oscillation along the y-axis (o s ) and the sliding motion along the boundary (o&), respectively. 
k 2 = a 2 characterizes the separatrix motion and the unstable isolated periodic orbit. An orbit on 
the separatrix alternately passes through one or the other focus between reflections. The unstable 
periodic orbit in the center of the separatrix performs an oscillation along the x-axis (o u ). It is 
the only orbit which passes through both foci between consecutive reflections. The length of this 
orbit is L u = 4, the stable orbit has L s = 4\/ 1 — a 2 and the length of the sliding orbit is just the 
circumference = 4£(a) of the billiard. £{k) is the complete elliptic integral of the second kind 
in the notation of |23|, 24]. The action S and the period T of these orbits can easily be obtained 
from S = V2EL and T = dS/dE = S/{2E) = L/V2E. 

The calculation of the action variables i" = (7 P ,/^) = § p P dp, ^ § P(pd(ft), the frequencies 
u = dH/dl, and the winding number w = co^/ujp of the elliptic billiard can be found in [10|. 
Since the system is invariant under reflections about the x- and y-axis, we will also consider the 
desymmetrized elliptic billiard, that is a quarter of the full ellipse. For type R (k 2 > a 2 ) the actions 
I and the winding number w = Cj^/ujp of the desymmetrized billiard are 

I p = ^ (sm X -K£( X ,Z)) = I p , 

k = = ±fo, (8) 




Figure 2: Energy surfaces (a) for the values E = 0.5, 1, 1.5, 2, 2.5, 3 and the winding number (b) of the desymmetrized 
elliptic billiard with parameter a — 0.7. The line k 2 = a 2 (I p — —^7^,) represents the separatrix motion and the 
unstable oscillation along the x-axis. The line k 2 — (1$ = 0) marks the stable oscillation along the j/-axis, and 
k 2 = 1 (Ip = 0) corresponds to the sliding motion along the billiard boundary. The dotted lines represent the real 
complex tori of type O c . 



with sin 2 % = (1 — k 2 )/(1 — a 2 ). For type O (k 2 < a 2 ) they are given by 



2 P ' 

\h , (9) 
w , 

with sin 2 ip = (1 — a 2 ) /(I — k 2 ). K{k) is the complete elliptic integral of the first kind, T{p,, k) and 
E(fj,,k) are incomplete elliptic integrals of first and second kind. 

Figure ^a shows the energy surfaces H(I) = E in action space. All these lines have the same 
shape, because the actions scale with \fE. We denote the curvature of the energy surfaces by C and 
its sign by (3. This means (5 = +1 for the patch being concave away from the origin (type R) and 
(3 = — 1 for the convex patch (types O). In the limit k 2 — > — oo the curvature tends to zero. The 
energy independent winding number w is shown in Fig. |2]b. The cusp at k 2 = a 2 has the limiting 
value «; = 1. The winding number of the stable isolated orbit o s with k 2 = is 2arccos (a)/vr. 

Beside the three special periodic orbits o s , o u and Of, there are families of periodic orbits on 
resonant tori. A /x-resonant torus is determined by its rational winding number w(k) = [i^j^p or 
equivalently by the frequency vector u being proportional to fi = (/i p ,/i^), where [i p and fi^ are 
relatively prime integers. We denote the action vector of a ^t-resonant torus by I 1 * . The action of 
a prime non-isolated orbit is given by S(I^) = Figure || shows prime orbits of (fi p , fj,^)- 

resonant tori for type R and O for the full elliptic billiard. According to (||), the winding number 
for type R is reduced by a factor 2 if compared with Fig. ||b. The orbits are chosen in such a way 
that they are always symmetric to the y-axis and, if possible, also to the x-axis. For type R the 
integer fi p counts the number of reflections at the boundary and fi^ the rotations about the origin. 
For type O the integer [i p gives half the number of reflections, and /i^ is half the number of times 
that the orbit touches the caustics cos^ = ±K/a. 

In the course of this paper it will be important to distinguish different kinds of non-real orbits 
which have to be taken into account in the uniformization. Since we are dealing with a separable 
system, it is sufficient to classify non-real orbits for one degree of freedom systems. We define 
and distinguish non-real orbits via the way a (generalized) action integral can be assigned to 
them. Allowing complex values for the momentum p and the position q, we are dealing with 
27r/ = <f p(q) dq, where p(q) 2 = 2(E — V(q)), and the integral is taken between (possibly complex) 
turning points, i.e. the (possibly complex) zeroes of E — V(q). A real orbit has an action integral 
connecting real turning points on a path (in the complex g-plane) for which p is real, i.e. just 
along the real g-axis. Of course the value of this integral is the same for any path in the complex 



w = W *)/£(*) 



R O 




Figure 3: Periodic orbit representatives of low (/i p , /^-resonant tori for a = 0.7, shown together with the corre- 
sponding caustic. For type 7? we have H^/Hp = w = w/2 < 1/2, while wn > 2arccosa for type O. 

plane encircling the two turning points. The special integration path chosen can be thought of 
as the orbit of the particle. A "tunneling orbit" also connects real turning points, however, the 
energy E is smaller than V(q), such that p is purely imaginary along the real (/-axis, hence giving 
a purely imaginary action. There are two more non-real orbits connecting complex turning points. 
The "scattering orbit" is related to quantum mechanical scattering resonances forming above a 
potential maximum. Its path of integration is taken to be an anti-Stokes line. It is the line of q 
in the complex plane for which the integral p{q') dq' is purely imaginary, where q\ is a turning 
point. In the case of an even potential with maximum at 0, as it is the case for the elliptic billiard, 
the anti-Stokes line is just the imaginary (/-axis, and the momentum is real. By construction the 
action of a scattering orbit is purely imaginary. In the following, tunneling and scattering orbits 
will be treated similarly, we refer to both of them as "imaginary complex orbits" . They are related 
to the "ghost orbits" described in |25|] and the barrier penetration integral of Miller [jl2[ gives i 
times their action (i = \/ — 1). The final type of non-real orbit is related to orbits below a potential 
minimum. They are obtained by integrating along a Stokes line (on which the integral is always 
real) between complex turning points, hence their action will be real. In the case of a symmetric 
potential with minimum at the Stokes line is just the imaginary (/-axis, and the momentum is 
purely imaginary. We call this kind of non-real orbit "real complex orbit", because its action is 
real. These orbits are the "complex orbits" described in [|l5| . 

In separable two degrees of freedom systems all the combinations of the four possibilities might 
occur in principle. Actually there is some arbitrariness in how to combine the possible real and 
non-real motions from the separated degrees of freedom. However, in the present case it seems to 
be natural to form pairs from tunneling and scattering orbits, giving rise to imaginary complex 
tori. These tori are denoted by TO and TR in Fig. [l|. Continuing the TO torus below the potential 
minimum gives rise to the TO c torus. The real complex orbit is the natural continuation of an 
elliptic orbit, which disappears at the minimum of the potential. The torus denoted by O c in Fig. [l] 
is therefore real in its p-part and real complex in its <?3-part. These real complex tori lead to an 
energy surface that extends beyond the positive quadrant, as discussed by Berry and Tabor. There 
are no continuations of the R and TR tori, since for the p- motion the hard billiard wall does not 
create any complex zeroes for low energy. Therefore the energy surface cannot be continued at 

K 2 = l. 

For type O c (k 2 < 0) the actions and winding number of the desymmetrized billiard are given 

by 

I p = r v / T^? sin r + v^^^r, A;)- £{r,k))) = \l p , 

k = ^ (v 7 ^^? (£(k) - £(*))) = , (10) 

w = J 7 (T,k)/fC(k) = w , 

with sin 2 r = 1— a 2 and k 2 = —K 2 /(a 2 — k 2 ). The formulas (|l0|) equal the formulas (|9|) for imaginary 



k. It is important to notice that although the actions are real, the tori are classically not allowed, 
because both position and momentum are complex. These formulas define the continuation of the 
energy surface shown in Fig. ^a. The winding number of the tori of type O c decreases very slowly 
to zero for k 2 — ► — oo. Resonant real complex tori can also be defined, because their real frequencies 
can fulfill a resonance conditions. The action S and the period of a non-isolated orbit on a resonant 
complex torus are always positive, although the angular component 1$ is negative. 



3 Exact quantum mechanics 



The Schrbdinger equation including the boundary condition separates in elliptic coordinates. Using 
the angle-radius parametrization (^, p) gives the standard form of the Mathieu equation 



d\G + (A - c 2 cos 2 4>)G = 
d 2 p F- (A-c 2 cosh 2 p)F = 0, 

which follows directly from (^||7|). The relation to the physical parameters is 

c 2 = 2Ea 2 /h 2 
A = 2Ek 2 /H 2 . 



(11) 
(12) 



(13) 
(14) 



For Dirichlet boundary conditions the eigenfunction ^(cf),p) = G((f))F(p) must be zero on the 
billiard boundary, which gives F(p max ) = 0. The solution in the angular variable must be periodic 
with period 2ir, to give a physical solution. Floquet theory guarantees the existence of solutions 



with period a multiple of 7r, because (11) is linear with 7r-periodic coefficients. In the classical 
terminology the special values of the parameters, for which ir or 2tt periodic solutions of (|TT|) exist 



are called characteristic values. For even solutions, i.e. G(0) = 0, they are denoted by a = A — c /2, 
and similarly by b for odd solutions. Since the Mathieu equation is an equation of Sturm-Liouville 
type, these eigenvalues are all real, ordered as ao < bi < a\ ■ ■ ■ for fixed c 2 and the corresponding 
eigenfunctions have i zeroes in the interval cj) E [0, 7r). Solutions with even i have period tt, those 
with odd i have period 2ir. 

If G(4>) is even, then ^ is symmetric with respect to the x-axis: ^(x, y) = ^(x, —y), see @. We 
denote this symmetry by ir y = +1, respectively by ir y = —1 for odd G with ^(x,y) = — ^>(x, — y). 
Similarly ^ can be even or odd with respect to x, which is denoted by ir x = ±1. The four possible 
parity combinations are given in Tab. |, together with the sign of the value of G(7r/2/c), k = 0, 1, 2, 3. 
In order to obtain a smooth wave function ^ in the part of the x-axis connecting the foci, the radial 
solution must satisfy -F(O) = if G is odd. If both parities are the same, the angular solution has 
period tt. In the last column we indicate which coordinate axis becomes a nodal line for the wave 
function of the corresponding parity. The given signs are only defined up to a global factor. 
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Table I: Symmetry properties of the wave function depending on the four parities. The four states which have the 
same quantum number are listed in the order of increasing energy. 



We define the radial quantum number r as the number of zeroes of F(p) in the range p G 
(0,Pmax)> i-e. not counting the zeroes at the boundaries. Similarly the angular quantum number I 
gives the number of zeroes of G(<p) in the range 4> E (0, 7r/2), again not counting possible zeroes at 



Figure 4: Contourplots of the probability density of the wave functions for the four different parity states (n x ,ir y ) 
for quantum numbers (r, V) = (0, 0) and with a — 0.7. The contourlines are equally spaced from to the maximum 
probability density of the wave function. 

the boundaries. This choice of quantum numbers is consistent with the EBK quantization (ffl) for 
the symmetry reduced system. Therefore it does not directly give the number of nodes of the wave 
functions in the full system: depending on the parity there are 0, 1, 2 or 3 additional zeroes in 

the range [0, tt] for the states with parity ++, — h, H — , and , respectively, such that the total 

number of angular nodes is 4/ + 2 — tt x — Tr y . The energy eigenvalues for the four states with the same 
quantum numbers are in the same order. In general we denote a state by (r, l) 7Tx7Ty . Note that the 
four different parity combinations correspond to the description of the symmetry reduced quarter 
billiard, where on the coordinate axes Dirichlet (parity —1) or Neumann (parity +1) boundary 
conditions are required. To illustrate the symmetries, the probability density of the (0, 0) state for 
each parity is shown in Fig. f|. 

The transformation cos 4> = r\ and cosh p = £ gives the algebraic form of the Mathieu equation 

(1 - r, 2 )d 2 g - rjdrjg + (A - c\ 2 )g = (15) 
(l-e)dlf-^f + (X-c 2 f)f = 0, (16) 

from which it is obvious that the angular and the radial equation are actually the same, only 
evaluated on different ranges of the independent variable, rj £ [—1,1], £ G [1, l/o]. This equation, 
although it has regular singular points at £,77 = ±1, is better suited for the numerical solution of 
the eigenvalue problem. The requirement for the solution to be smooth at these points replaces the 
periodic boundary condition. 

The Mathieu equation is a special case of the spheroidal wave equation, 

(1 - x 2 )d 2 h m - 2(m + l)xd x h m + (A - 1/4 - m(m + 1) - c 2 x 2 )h m = , (17) 

which appears in the case of the billiard inside the rotational symmetric ellipsoid; m is the quantum 
number of the angular momentum of rotation around the axis of symmetry. With m = —1/2 in 
( |l~7|) we reobtain (|l^) and ([l~6|). We will see that m = +1/2 also produces solutions of the Mathieu 
equation. 

In order to obtain solutions with 7r y = —1, which must be zero at </> = 0, tt, we must construct 
solutions which are zero at the regular singular points. It is necessary to factor out this behavior 
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Figure 5: The solutions of (jl^) for the state r = 1, I = 4 with the four possible parities, i.e. even or odd at x = 
(iv x ) and m — —1/2 for ir y = 1 or m = 1/2 for ir y — —1; a = 0.7. 

by the ansatz g(x) = yY— x 2 /i!/ 2 (x), /(x) = \/x 2 — 1/11/2(2:), which transforms ([H]) and ([!(]) into 
the equation 

(1 - x 2 )8 2 x h 1/2 - 3xd x h 1/2 + (A - 1 - c 2 x 2 )h l/2 = 0, (18) 

which again is a special case of the spheroidal wave equation (|l~7|), in this case corresponding to 
m = +1/2, however. 

In Fig. H four solutions of ( JT7| ) are shown. The radial and the angular part smoothly join at 
the regular singular point at x = 1. We conclude that the spectrum of the 2-dimensional billiard 
is obtained from the spectrum of the 3-dimensional rotational symmetric billiard if the angular 
quantum number m in the spheroidal wave function is set to a half integer "spin" number ±1/2 
instead of to an integer number as in the 3-dimensional problem. 

In the standard theory of the Mathieu equation c is fixed and the eigenvalue A is determined. 
In the billiard problem we have to simultaneously satisfy also the boundary condition for the radial 
equation. Since both separation constants A and c 2 appear in each equation, although the variables 
are separated, the separation constants are not separated, see e.g. p^| , which requires a nonstandard 
approach to the numerical solution of this Sturm-Liouville eigenvalue problem. For the angular 
equation (|l5|) we have boundary conditions at r\ = —1 and at r\ = 0, corresponding to <fi = tt and 
(j) = 3tt/2, for the radial equation at £ = 1 and £ = 1/a, corresponding to p = and to the billiard 
boundary. Introducing a new independent variable by £ = 1/a + C(V a — 1)> t ne ranges of £ and 
j] coincide, and we can use a standard shooting method as, e.g., described in pTj] , to solve both 
coupled eigenvalue problems simultaneously. Even though the final solution is smooth, as shown in 
Fig. [|, it is not possible to numerically integrate through the singularity. Instead the integration 
always starts an epsilon away from the singularity with an analytically calculated regular initial 
velocity In order to find a specific state it is necessary to have a fairly good initial guess for 
the eigenvalues, otherwise the shooting method might converge to another state. We calculate the 
initial guess via the uniform WKB approach described in the next section, and have found that 
this always works. The results of these computations are shown in Fig. |(| and some eigenvalues 



are listed in Tab. Ill 



In Fig. g a few symmetric wavefunctions are illustrated. The wave functions for the other 
parities are related to the ++ states show in Fig. R as illustrated in Fig. El. In the left column the 
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Figure 6: The eigenvalues (E, k 2 ) of the eigenstates of the elliptic billiard with a — 0.7. The region of type R states 
(rotational, k 2 > a 2 ) is shown in light grey, type O states (oscillating, k 2 < a 2 ) in dark grey. 
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Figure 7: The classical actions J = 27 in units of h corresponding to the exact quantum mechanical eigenstates in 
the elliptic billiard for a = 0.7. Except for states close to the separatrix, they are located on a lattice given by EBK 
quantization. The structure near the separatrix can be explained by uniform WKB quantization. Grey code as in 
Fig. |. 



Figure 8: Contourplots of the probability densities of the wavefunctions (r, l)++ of the elliptic billiard with a = 0.7. 
The vertical axis shows r = 0, . . . , 5 and the horizontal axis shows I = 0, ... ,5. The classical caustic corresponding 
to the second eigenvalue is also shown. 

localization around the stable isolated periodic orbits is clearly visible. It becomes stronger when 
r is increased. In the bottom row the same happens for the orbit o^. Both cases can be explained 
by considering the WKB wavefunctions for the corresponding states, which is exponentially small 
outside the classical caustics and nonzero inside. Although the WKB wavefunction diverges at the 
classical caustic it is correct insofar as the quantum probability density is relatively high close to 
(and inside) the classical caustic, as can be seen e.g. for the (4,4)++ state. The surprising fact, 
which can not be explained by the above reasoning, is the localization around the focus points, 
e.g. in the (1,4)++ state. This phenomenon occurs at the transition from the rotational states 
(classical caustic type R) to oscillating states (classical caustic type O), i.e. close to the classical 
separatrix. Eigenstates with higher quantum numbers can show even stronger localization around 
the foci, because their second eigenvalue can be found closer to the critical value k = a = 0.7. 
Besides this very strong localization on the focus points these states also show a "scar" of the 
unstable orbit, e.g. see the (2,5)++ state. Although the probability density is low compared to 
the one in the focus points, it is high compared to the region further away from the x-axis. The 
behaviour outside the foci can again be explained by the WKB wave function, which diverges for 
y = for the critical value k = 0.7. We conclude that in the integrable elliptic billiard states 
localize around stable periodic orbits, and they also show "scars" along the unstable periodic orbit, 
with an additional strong localization on the focus points of the ellipse. 

The oscillating states of type O are always non-degenerate. The rotational states of type R 
become more and more degenerate when the distance from the separatrix is increased. The classical 
reason for this increasing degeneracy is the fact that there are two tori for fixed constants of motion, 
connected via time reversal symmetry. In a simple EBK quantization these states would be exactly 
degenerate. Two type R states with the same r and I become (approximately) degenerate if they 
have period 2ir, i.e. if TT x TT y = —1. If, however, the period is ir, ir x Tr y = 1, two states with the same 
r but I differing by 1 become degenerate. This is consistent with the fact that the total number of 
nodes on the circle 4> £ [0, 2tt) is 4/ + 2 — ir x — ir y , which gives the same number for the degenerate 



states. The shift from degenerate states of type R to non-degenerate states of type O can most 
clearly be seen if the exact quantum eigenvalues E, k 2 are transformed into the classical action 
space using equations (g) and @, see Fig. [?]. In order to represent the actions of the full system 
for both types of motion in one picture in such a way that on average one state occupies a box 
of size H 2 , we use the doubled action of the symmetry reduced billiard J = 21. This amounts to 
doubling the action I p for type R tori, which can be understood from two facts: firstly there are 
two classical tori with the same action for type R tori and secondly the corresponding quantum 
mechanical states are degenerate in the EBK approximation. The type of movement that each 
parity state can perform in its "semiclassical quantum cell" will become clear in the next section. 



4 Uniform WKB quantization 

Consider a Hamiltonian H = p 2 /2 + V(q) with a potential as e.g. given by the effective potentials 
of sec. |2|, see Fig. |l[ Let V(q) < E for q < q\ and q > q^ and V(q) > E in the range q\ < q < q%. 
The WKB solutions to the left (j = 1) and to the right (j = 2) of the potential barrier are 

^(q) = (A+ expiiSjiqyh) + Aj expHS^/ft)) /VRq), (19) 

where 

p(q) = s/2(E-V(q)), Sj(q) = f p{q') dq' . (20) 

J qj 

The matrix connecting the constants Af and A^ to A^ and A^ (see e.g. the review of Berry and 



Mount [14] and the references therein) is given by 



f A+\ ( A+ \ n/h ( Vl + e~ 2 ®/ h -i \ 

Ur M U)' I * vTT^J' (21) 

with the penetration integral = —ij^p(q')dq'. In the terminology of sec. || this is i times the 
action of an imaginary complex tunneling orbit. This formula is also valid in the case where the 
energy E lies everywhere above the potential V(q). Then the classical turning points become com- 
plex (qi complex conjugate to 52) and the penetration integral 6 becomes negative, corresponding 
to an imaginary complex scattering orbit. In a potential with several turning points in each region 
away from the classical turning points a WKB wave function is reasonable. In classically forbidden 
regions the phase S becomes complex leading to a real exponential. Uniform semiclassical quantiza- 
tion conditions may be obtained by a piecewise connection of these wave functions, and by imposing 
the correct boundary conditions, e.g. an exponential decay in classically forbidden regions, periodic 
boundary conditions in the case of a rotor with a potential, or Dirichlet conditions in the case of a 
hard potential wall, e.g. in a billiard. The resulting quantization conditions take into account real 
tori and imaginary complex orbits. 

The quantization conditions for the planar elliptic billiard are obtained from the effective po- 
tentials shown in Fig. [l]. We allow for negative values of p, which can formally be achieved by 
letting <p vary between and ir only. This is helpful in order to incorporate the different parities 
(see the boundary conditions in Tab. ||) in the semiclassical wave functions Fj(p) (Eq. (|l9| ) for the 
radial degree of freedom). The wave functions Fj(p) are symmetric or antisymmetric with respect 
to p = 0, corresponding to 7r y = +1 or 7r y = — 1 (see Tab. |), consequently 



At 
A~ 2 




(22) 



The matrix P describes the phase shift along the classically allowed region 




where it is important not to use the action I p of the full billiard, because this quantity jumps by a 
factor of 2 as the effective energy increases beyond the top of the barrier. Instead we use J = 21 
as introduced in the last section, such that the phase shift is continuous at the transition from type 
O to type R motion. Inserting the Dirichlet boundary conditions -Fi(po) = ^2(^3) = into (|l9|) 
leads to 




oc P 



1 



1 



a; 



oc P- 



1 



■1 



(24) 



The fourth relation among the constants Aj is given by (|2l|), where the barrier penetration 
integral = Q p is calculated from one turning point to the other and is given by 



-2V2Ek (£{^l-a 2 /K 2 ) - K{^/l-a 2 /K 2 ) 
-2V2Ea (£(y/l - K 2 /a 2 ) - ^K(^l-K 2 /a 2 
-2\/2Eyf^ 2 



K 2 £(a/ya? 



type TR 
type TO 
type TO c 



(25) 



The action of imaginary complex orbits of type TO c will be needed in the next section. 

Composing the connection formula (pi]), the symmetry (p2[ ) and the boundary conditions (p4|) 
gives 

' 1 \ f -l\ 



PUP 



1 



1 



(26) 



Decomposing this complex Equation into its real and imaginary part we obtain the quantization 
conditions for the radial degree of freedom, 



and 



cos (nJp/h) 



sin (irJp/h) 



-7T„ 



(27) 



(28) 



y/l + e -2S P /h 

which have to be fulfilled simultaneously. The second Equation is not independent of the first one, 
it just selects half the number of the solutions of the first Equation. 

For the WKB wave function G(<j>) we use V e s((j)) (Fig. |l|b), and we impose periodic bound- 
ary conditions. Following the calculations of Miller |l2| for the different parities we obtain the 
quantization conditions 



COS 



(29) 



and 



sin (ttJ^/H) 



7T, 



+ e 



(30) 

is a result of the especially 
just selects the x-parity 



with = — Q p . The equivalence of the absolute values of P and 9 
symmetric separation of the Hamiltonian (|5|). The condition in Eq. 
of the solutions of fl29|). 

In Table || the limiting cases for large |0| of the right hand sides of the Equations (27-3C) 
are given. The resulting EBK quantization conditions for the full billiard with quantum numbers 
(rip, riffr) and the connection to the quantum numbers (r, I) of the symmetry reduced billiard are 
indicated. For type R motion these conditions are equivalent to the EBK quantization conditions 
for the circular billiard. 

How to calculate the (E,k 2 )- or the ( J p , J^-spectrum from the conditions in Eq. (p?|), (pS|), 
(|29| ) and (|30|)? The actions and the barrier penetration integrals are functions of (E,k 2 ). With 
Newton's method we can numerically obtain (E, k 2 ) as functions of (J p , J$) and therefore P and 
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Table II: The limiting cases of the right hand sides of the Equations ( p7| - |3(i| ). The quantum numbers for the 
corresponding EBK quantization conditions are n p ,n^ £ NU {0} for the full system and r,l S N U {0} for the 
dcsymmetrized billiard. 



00 as functions of (J p , J^). Ignoring (pq) and (|30|) for the moment, the remaining Equations (p7|) 
and ( |29| ) can thus be rewritten in the form 

/^(Jp,^) := cos (n J p / h) +TT y g 1 (Q p (Jp,J^)) = 
f2 X,7Ty (Jp,Jcf>) ■= cos (nJ^/h) - ir x Tr y g 2 (@ ( i > (Jp, J,p)) = 0, 



where gi and 52 are smooth functions onto the intervall (0, 1). A solution of ( |3l| ) can be viewed as an 
intersection of the lines y = and f^ x,7Ty = in action space with the same corresponding parities. 
The solutions lie inside a box with the edges defined by the extremal values of g\ and g 2 - Because 
of the periodicity of the left hand sides of Equations (27-|30|) these boxes are always arranged in the 



same way inside a cell in the ( J p , J^-space. We call these quadratic cells "semiclassical quantum 
cells" with quantum numbers (r, I). Their edges have length AJ p = AJ^ = 2h and they tessellate 
the whole (J p , J^)-space. Inside a quantum cell there are four quantum states, one for each parity. 
Each state of fixed parity is confined to a "parity box" of width h/2, shown by bold lines in Fig. |9[ 
The size of the parity box is a result of the fact that the Maslov indices change by two upon the 
transition of the separatrix, see Tab. |I| These parity boxes take care of the remaining conditions 
(|28| ) and (|30|): inside a box they are automatically fulfilled when (^) holds. To find the solution 
guaranteed to lie in the interior of a box, we use the bisection method described in |p9[ |. Having 
got all solutions in action space, we compute again the solutions in (E, K 2 )-space with Newton's 
method to obtain the semiclassical eigenvalue spectrum. Introducing the semiclassical quantum 
cells in action space ensures that always the right state is found, even when two states are almost 
degenerate. The method can easily be extended to the case of more than two degrees of freedom. 



Figure 11 and Tab. Ill show the result of the semiclassical calculation compared to the exact 



quantum mechanical results. The distinction between the parities is omitted in Fig. 11 because this 



becomes already clear from Fig. ^. In Fig. 1C the relative errors of the semiclassical eigenvalues 
(-Eqm — E sc )/E qra and (K qm — K sc )/K qra are plotted versus K qm , such that the behavior of the error 
with respect to the position of the state on the classical energy surface can be seen. The semiclassical 



energy eigenvalues are almost always too low, see also Tab. [II. The opposite holds for the values 
of the second eigenvalue, which usually are too high. The only exception (for both eigenvalues) 
occurs in the neighborhood of the classical separatrix. The fact that Fig. |l0| looks rather symmetric 
indicates that the relative errors of the two eigenvalues are strongly correlated. Concerning the 
semiclassical limit we see series of states with increasing r quantum number and decreasing error 
for small ft, and series with increasing / quantum number and decreasing error for large k. However, 
increasing / with k, < a does not decrease the error. The error of the whispering gallery states (0, 1) 
is particularly large. In general, however, the agreement between the exact eigenvalues and the 
semiclassical values is very good. 

Now we study the change of the regular lattice obtained from EBK quantization induced by the 
uniformization. There are three typical situations in action space represented in Fig. region R 
(a) and region O (b) both far away from the separatrix and the region close to the separatrix (c). 



In region R the EBK states (r, / + 1)++, (r, /) (e.g. Tab. Ill rows 9 and 10) and (r, l)^ , (r, /) 
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Figure 9: Schematic view of the four states in a single semiclassical quantum cell (AJ P = AJ^ = 2h) for (a) the 
region R, (b) the region O and (c) close to the separatrix corresponding to the states with quantum numbers (0, 1). 
The intersections of the lines f* y = and /^'"^ = with equal parities (dots) are the solutions of the quantization 
conditions (|27|), (pg|), ( p9| ) and ( pSC| ) if they are located inside a parity box. The squares and the triangles are the 
states obtained from EBK quantization in regions R and O, respectively. Parity boxes that are crossed by the bold 
dashed separatrix may contain two EBK states (bold squares and triangles) or no EBK states, because both corners 
(the thin squares and triangles) are in the wrong region. The EBK Maslov index is (3, 0) in region R (where J p = 2I P ) 
and (4, 2) in region O. 




Figure 10: The relative error of semiclassical energy eigenvalue (circles) and k eigenvalue (pluses) plotted versus 
k. The obvious patterns result from series of states with the indicated quantum numbers. The arrows point into the 
direction of increasing quantum number. 



(e.g. rows 16 and 17) are almost degenerate. This corresponds to the degeneration of clockwise and 
counterclockwise rotations in the planar circular billiard. The EBK states form regular lattices, with 
mesh size A J p = 2h and AJ^ = h with two states on each corner in region R and AJ p = A = H 
in region O. Far away from the separatrix the uniform semiclassical states lie near the EBK states 
with Maslov index (3, 0) for type R and (4, 2) for type O. Approaching the separatrix a transition 
between the two regular EBK lattices takes place. Crossing the separatrix from region R to region 
O the shift of the EBK states is 6J P = —H/2 and 5J^ = h/2 for ir y = +1 and 5J P = h/2 and 
8J^ = —h/2 for ir y = —1. The uniform quantization smoothly joins the EBK lattices along the 
separatrix. It is natural to interpret the location of the semiclassical states as two separate meshes: 
one for ix y = +1 and one for ir y = —1, indicated by full and dotted lines in Fig. 11. 

Figure |9|c shows that close to the separatrix not only the accuracy of the EBK states is unsat- 
isfactory but also that EBK quantization sometimes yields a wrong number of states. In Fig. [9|c 
the situation for the (0, l)-quantum cell is shown. The EBK states are marked by bold squares in 
region R and by bold triangles in region O. A square represents two EBK states because of the 
degeneracy in region R. The +— parity box contains one EBK state, marked by the bold triangle. 
The opposite corner of this parity box could be reached with the Maslov index defined for region 
R. However, it is located on the wrong side of the separatrix, and therefore not a true EBK state 
indicated by a thin square. By the same reasoning the parity boxes ++ and — + contain no true 

EBK state. On the contrary the parity box contains two true EBK states. This indicates the 

main deficiency of EBK quantization for systems that exhibit a separatrix. 




Figure 11: Comparison of the exact states (circles) to the semiclassical states (intersections of the lines f^ v = 
and f 2 x ' y = with the same corresponding parities). The dashed line is the energy surface E — 200. J p and J$ are 
measured in units of h. 



5 The Berry- Tabor trace formula 

In the last section we have seen how the imaginary complex orbits lead to a Maslov index smoothly 
varying across the separatrix. The real complex orbits can be taken into account by a completely 



different approach introduced by Berry and Tabor [15]. The goal of this section is to incorporate 
the imaginary complex orbits into the Berry- Tabor trace formula, such that all kinds of classically 
forbidden tori are taken into account. 

Starting from EBK quantization Berry and Tabor derived a formula for the semiclassical density 
of states 

n sc {E) = ^{E-H({n + j)h)) (32) 
n 

in terms of resonant tori for an integrable system with / degrees of freedom. They transformed 
Eq. (32) via the Poisson summation formula and performed the integrals in stationary-phase ap- 



proximation. The result is 

2 ^ I ^ cos (q(S(I»)/h - § an) + f /?) 



'Is 



(E) = n(E) H j^j- > /-i . 2^ — • ( 33 ) 



9=1 Q 



The first summation is over all relatively prime non-negative integers fi, i.e. over families of 
prime orbits, while the second summation runs over all their repetitions. The Maslov indices are 
always the same since the energy surface is assumed to have no separatrix. The term /x = is 
excluded from the sum and denoted by n(E); it gives the mean density of states, the so-called 
Thomas-Fermi term. For the elliptic billiard we have n(E) = Ab/(27rh 2 ) according to Weyl's law 
where = 7r\/l — a 2 gives the billiard's area. In [^] it is shown that next order correction to 
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Table III: The quantum mechanical eigenvalues (-Eqm,Kq m ) and the semiclassical eigenvalues (E ac ,K^ c ) of the 
billiard in the ellipse for a = 0.7 with E qln < 26 or 1000 < £ qm < 1030. 



Weyl's law, proportional to the circumference of the billiard, is contained in the whispering gallery 
orbits in the remaining sum. 



To improve the convergence of the series (33) it is advantageous to introduce a smoothed density 
of states n sc ([E; 7) by giving E a small imaginary part £7. Then a <5-peak changes into a Lorentz 
function with half width energy 7. The semiclassical formula for n sc (E; 7) differs from formula ([33|) 
by a decay factor D = exp (— r yqT(I fl )/h) before the cosine term in the second summation. 

Berry and Tabor show that formula ( |33|) can be improved by taking into account resonant 
real complex tori. This removes unphysical discontinuities in the density of states resulting from 
contributions of resonant tori which are suddenly classically realized as E changes. For a system 
that scales with respect to the energy this cannot occur, nevertheless the system parameter a can 
take over the role of the energy, e.g. for the elliptic billiard all real complex orbits become realized 
for a — ► 1. Even without a parameter these corrections are sensible because they improve the 
situation when the stationary phase approximation is bad because there is a stationary point close 
to, but outside the range of integration. The resulting trace formula incorporating real complex 
orbits for two degrees of freedom systems reads 

n sc (E; 7) = n(E) + JhE — fP= £ D 

V 1 y\ / sin (q(S 2 a fl)) Bin(g(Si/fi-fa/J)) \ 

4^ \i.KT»\\J\,,\\ntT^\\ ^ I A 2 ^ h x Jq J (34) 
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,_ J_ yy y. ( sin(g(5 2 / h- § a 2 (J,)) _ sin(q(S 1 /h~§a 1 fJ,)) \ 



M 9=1 

where £1,2 are the actions, 01,2 the Maslov indices and i\ 2 the normalized derivatives with respect 
to a second constant of motion £, which parametrizes the energy surface in action space; all quan- 
tities with indices 1,2 are evaluated at the boundaries £1 and £2 of the energy surface. A is the 



amplitude and 6 the argument of the complex Fresnel integral 



' dxe^ 1 , (35) 



2-7T JKi/Vh 

where A 1)2 € R are defined by A? 2 = 2q/f3(S 1>2 - S) and Ai >2 < if £ lj2 < £ M , where = 1(f). 
The series ( |34|) holds for real and real complex tori. 

Formulas (33) and fl34| ) can be applied to simple systems [3(J but not to systems with 



separatrices like the elliptic billiard. In the derivation of Berry and Tabor the energy surface is 
assumed to be strictly convex or strictly concave and to have a smooth curvature. This is not the 
case for a system with separatrices, e.g. the elliptic billiard, as can be seen in Fig. 0. We are going 
to study three versions of the Berry- Tabor trace formula with sucessive improvements concerning 
states close to the separatrix. A first approach, similar to [^lj, would be to divide the energy 
surface into smooth patches and to consider (|3^) or ( |34| ) for each patch. We call this the EBKBT 
quantization. The manifestation of the separatrix in the EBK approximation as a discontinuous 
quantization condition is carried over to the EBKBT quantization by the appearance of the wrong 



number of terms in the summation (33) or (|34j). In the semiuniform WKBBT quantization we 



correct this deficiency by taking the Maslov index in (32) as a uniform Maslov phase, smoothly 
varying across the separatrix. However, the new term is considered constant in the stationary 
phase approximation. This improves the results and still gives a sum over periodic orbits. The 
best results are obtained if the Maslov phase is fully taken into account in the stationary phase 
approximation, leading to summation over (in general) non-periodic orbits in a uniform WKBBT 
trace formula. 

The resonant tori of an integrable two degrees of freedom system are characterized by the 
winding number w{k^) = \x^l 1 [i p . For a billiard system without potential w is independent of the 
energy. Thus we can obtain all resonant tori for one reference energy. The actions, the frequencies 
and the reciprocal curvature scale with \[E. Hence all these quantities must only be calculated 
for the reference energy. For the elliptic billiard the winding number w is restricted to the interval 
(0, 1). Thus for resonances up to a fixed order we have < \i p < /x max and < fi^ < fj, p . For 
the summation over q, i.e. over the repetitions, we incorporate all orbits with period qT less than 
a cutoff time T c = ^ In | to guarantee that the error in the second summation of formula (j33|) 



7 

is smaller than e. In the following we choose a = 0.7, h = 1, 7 = 0.3 and e = 10~ 8 . Thus the 
cutoff time is T c ~ 63.7. For numerical reasons it is not possible to approach the cusp of w (see 
Fig. ||b for w) arbitrarily close. We can only include all resonant tori with \k 2 — a 2 \ > 10 -14 , 
i.e. w < 0.472 for type R and w < 0.943 for type O. But the contribution of tori with larger 
winding numbers are very small because of the divergence of the curvature at k 2 = a 2 . Taking 
Mmax = 250, we include 26 395 and neglect 2 026 resonant tori in the summation. The remaining 
interval [0, 2 arccos (a) /it) ~ [0,0.506) carries the real complex tori (type O c ). Since for a = 0.7 
the real complex torus with w = 1/2 has the predominant effect we found it sufficient to consider 
resonant real complex tori with k 2 > K^ in = —15, i.e. w(k) > w(K m i n ) = 0.285. Resonant tori 
with k 2 < fv^in give negligible constributions, even though the curvature is very small, due to the 
destructive interference in the Fresnel integral (|35|). The number of incorporated real complex tori 
is 4205, 5 1 are left out. 

In the EBKBT approach we set a = (3, 0) and = 1 for the patch R and a = (4, 2) and 
(3 = —1 for patch O. On patch R the EBK states are exactly degenerate, thus all the contributions 
are multiplied by 2. Since there are only real complex tori of type O c , in Eq. ( |34|) for patch O we 
only consider terms corresponding to the boundary £1 (k 2 = 0) and no boundary terms for patch R. 



Figure 12 presents the spectrum n sc {E;^) calculated with the EBKB T quantization in the range 
12.5 < E < 47.5. For comparison the exact spectrum smoothed in a Lorentzian manner is shown. 
We focus on the states shown in the quantum cell (0, 1) close to the separatrix shown in Fig. U. The 



other states in Fig. 12 are reproduced fairly well. First we observe that the (0, 1)_| state is also 



reproduced quite well, which is due to the fact that the corresponding parity cell correctly contains 




Figure 12: Comparison of the exact quantum mechanical density of states (dotted line) to the semiclassical density 
n sc (_E;7) calculated from EBKBT quantization (solid line). The labels refer to the exact peaks. EBKBT peaks 
connected to parity boxes with the wrong number of states (see Fig. 0) do not give satisfactory results. 



one EBK state only. The peak for the parity box (0, 1)++ splits into two peaks with approximately 
half amplitude. The reason is that there is no true EBK state in this parity box. However the two 
corners of the box carry EBK states located in the wrong region (the thin square and triangle in 
Fig. ^) , and they still give a small contribution in the trace formula. The situation is similar for the 

(0, 1) y state. Here the EBK peak also splits into two small peaks; one accidently overlaps with 

the (1, 0)++ EBK state and the other one is seen as a small bump on the left side. In both cases the 

EBKBT quantization almost fails to produce a state. The parity box (0, 1) in Fig. || contains 

two EBK states. Since both of them are in the correct region (the bold square and triangle) they 

both produce peaks with almost the usual amplitude, which can be clearly seen by considering 

parity only, i.e. by looking at the billiard in the quarter of an ellipse. In Fig. 12 we instead observe 
that the (0, 2) ++ peak has almost doubled amplitude: this is the result of the assumed degeneracy 
of states in the region R. So here the EBKBT quantization produces one state too much. The 
corresponding peaks in the spectrum are called "large spurious peaks" in |^jf |. 

To remove these problems the tunneling through the potential barriers should be incorporated 
by a summation over imaginary complex orbits. For one-degree-of-freedom-systems the consider- 
ation of tunneling orbits can be done in a simple way |32|, 33 1. For systems with more degrees 
of freedom the problem is rather involved. In particular a Green's function approach in terms of 
action angle variables is still out of reach except for special systems |34|| . Thus we again start 
with Eq. (^) as Berry and Tabor did. At the heart of their approach are the EBK quantization 
conditions of the form (||) which allow for a resummation of the semiclassical density of states via 
Poisson summation. At this stage we allow for a varying a in Eq. (|l|), i.e. we base our calculations 
on a uniform quantization condition. Note that this uniformization is completely different from the 
uniformization of Berry and Tabor, which involves only real complex tori. We incorporate both 
approaches, with the result that all kinds of classically forbidden tori described in sec. [2] are taken 
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Figure 13: Comparison of the exact quantum mechanical density of states (dotted line) to the semiclassical density 
n sc (i?;7) calculated with semiuniform WKBBT quantization (solid line), semiuniform WKBBT peaks corresponding 
to parity boxes with no EBK state at all ((0, 1)++ and (0, 1)- + ) are improved as compared to Fig. [l2| The parity 
box (0, 1) with two EBK states still produces two peaks. 



into account. 

An effective Maslov phase a (I; tt x , Tr y ) varying smoothly on the energy surface can be extracted 
from Equations (27-|30|). Inserting the EBK quantization condition (|l]) into these Equations the 



quantum number n drops out because it appears as Inn in the argument of the trigonometric 
functions. Thus we obtain the effective Maslov phases as 

a p (I; tt x , n y ) = 7t y - arctan ( e epiI '" I ^ /h ) + 3 - -n y (36) 

and 

a<f,(I) tt x , TT y ) = TTy arctan ( e e 4V«>)M) + 2 - % - tt x . (37) 
Then we introduce new variables, behaving like continuous quantum numbers 

n = I/-Ia(J) (38) 

for each parity. The energy surfaces in these variables (one for each parity), which now depend 
on the energy through a, are shown in Fig. 14. The energy surfaces with equal n y have the same 



shape, but differ in a horizontal shift by the constant 1/2. For h — > a becomes a step function 
such that except for a shift and the parity splitting the classical energy surfaces are reobtained in 
the semiclassical limit. The quantization conditions for these variables read h £ N 2 and are trivial. 
In sec. §we studied how the lattice of states changes across the separatrix, while the energy surface 
was the same for all energies. Now we turn the point of view and introduce new variables h giving 
a trivial lattice but a more complicated energy surface instead. The lines of constant h p , in the 



usual action variables are shown in Fig. 11 



Since the quantization conditions in the new variables are of EBK type, we can repeat the 
derivation of the trace formula by Berry and Tabor. The phase being approximated in this deriva- 
tion is fin. Without separatrix the only non constant part in n is /, and the stationary phase 




Figure 14: The energy surfaces in the variables n for the energies 25, 300, 1026.76. In contrast to Fig. [H] here the 
shape of the energy surfaces depends on the energy but the semiclassical states form a regular lattice. For the energy 
E = 1026.76 the almost degeneracy of the semiclassical states (0, 16) _| and (0, 16) \- is shown. 



condition leads to resonant tori. In order to keep the summation over resonant tori we must assume 
that a(I;TT x ,TT y ) varies slowly in action space, which is also the approach followed in IS]. This 
gives a trace formula with two minor modifications: Firstly each sum is subdivided into separate 
sums for each parity. Secondly a now is different for every resonant torus. We call this method 
semiuniform WKBBT quantization. Figure ^ shows the semiuniform WKBBT spectrum for the 

same parameters as above. The states (0, 1)++ and (0, 1) y. (parity cells without EBK state) are 

reconstructed in a satisfactory manner. The (0, 1) state (parity cell with two EBK states) has 

improved, but it still contains two contributions, which is due to the fact that we assume a to 
be constant in the stationary phase approximation, which is a particularly bad assumption in the 
neighborhood of the separatrix. 

Thus we are finally led to the necessity to fully take into account the variation of a. In essence 
this means to take the surfaces in Fig. 14 as new energy surfaces, and define all quantities (most 
notably the winding number) with respect to them. Most important, the stationary points are now 
given by 

dl h da. & 

^ = 4^ = ^she ( ^-^- (39) 

Although far away from the separatrix |0| is large and this condition almost reduces to the ordinary 
resonance condition, in general Equations (^) and ( |34"[) are no longer summations over resonant 
tori and consequently the resulting semiclassical density of states is not determined by periodic 
orbits. Instead the tori to be taken into account in the sum are given by the rational values of a 
new effective winding number 

<«» 

such that the solutions of Eq. (|39| ) are given by w n (E,K) = \i^l\x p . The overall structure of the 
trace formula does not change, however, the frequencies and the curvature have to be replaced by 
the respective expressions obtained from the surface n(«). An important difference between the 
winding numbers w and w n is that the latter depends on the energy. This leads to an enormous 
numerical effort for the calculation of n sc (E; 7), because it is necessary to determine the stationary 
points for every energy separately. We call this method uniform WKBBT quantization in contrast 
to the semiuniform WKBBT quantization, since it neglects the varying a in the stationary phase 
approximation. The uniform WKBBT spectrum calculated with the same parameters as above is 
shown in Fig. 15. As expected the splitting of the (0, 1)-- state is now removed. In the evaluation 
of the sum we define the boundaries of the energy surfaces as = 0) and n(n = 1). This choice 




Figure 15: Comparison of the semiuniform WKBB T density of states (dotted line) to the density calculated with 
the uniform WKBBT quantization (solid line). The density obtained from the energy eigenvalues calculated in sec. ^| 
(dashed) gives the expected location of the peak, which is slightly displaced from the peak in the exact quantum 
mechanical density (long dashed). The main point is that the splitting of the (0, 1) — peak is removed. 



is somewhat arbitrary since one could also define the boundaries as h p = and = 0. But we 
want to focus onto states close to the separatrix where this arbitrariness is not relevant. The little 



peaks in Fig. 15 are artefacts which could be removed by including larger /x in the sum. 
The derivation of the effective Maslov phase ( p6| ) and ( p7| ) makes it obvious that in the uniform a 
the effects of tunneling and scattering orbits are incorporated into the Berry- Tabor sum. Now we 
want to show that the consideration of an effective a is equivalent to an explicit incorporation of 
imaginary complex orbits with constant a. 

We consider two crossed double well potentials, which is slightly more general than the elliptic 
billiard since the barrier penetration integrals 0i and 02 are then unrelated. The effective & 
obtained by the method described in sec. [| can always be written as <5j = —2Q/t: + a c %, with 
Qi = arctanexp (— 0j) and a c = (&&, a C 2) is a constant integer vector. We rewrite the cosine term 
in Eq. ( |33"1) with x = q{S/h — 7r/2a c fi) + 7r/4/3 and % = qfi{ > as the real part of 

e *(*+«Ci+92C2) = e ix( ri + itl yi ( r2 + it2 y2 ( 41 ) 

with ^ ^ 

r; = cos Ci = = and ti = sin 6 = = . (42) 

Vl + e- 2e * Vl + e 2e * 

Now rj and ti can be interpreted as the absolute values of the refection and transmission coefficients 
of the potential barriers for one reflection or transmission (|35|, [3(|), respectively. Expanding the 
last term of Eq. ( |4l| ) gives 

91 / \ 12 



k 1 ,! 1 =1 fc 2 . i 2 = 1 



V ^^^i Y: (Trrh l i, (43) 



where the summations go above all permutations of r±, t\ and r2, t2, i.e. we suppress the binomial 



coefficients. Finally we expand Eq. (^) and take the real part to give 



£ ©©-^i^y****- (44) 

fcl+!l=91 

Again the summation is above all permutations of n, ii and r2, £2- The phase in the cosine 



term in Eq. (44) can be interpreted as accumulation of phase shifts of tt/2 resulting from single 
transmissions. Thus Eq. (B) tell us how to extend the summation (^) to a summation over 
imaginary complex tori. One must incorporate all combinations of real and imaginary complex 
orbits with the same resulting period by replacing the cosine term in Eq. ([33]) by Eq. (f44|). 

6 The length spectrum 

Instead of calculating the density of states by a summation over resonant tori one can turn the 
tables and look at the length spectrum, i.e. the Fourier transform of the oscillating part of the 
density of states n osc (E) = n{E) —n(E). In the context of hyperbolic systems this viewpoint is 



referred to as inverse quantum chaology, see, e.g. [17, 18]. 

The familiar way to discover the appearance of periodic orbits in the quantum mechanical 
spectrum is to calculate its power spectrum. Similar to the phases in the Gutzwiller trace formula, 
the phases in the Berry- Tabor summation over resonant tori are proportional to the action of 
periodic orbit representatives of the tori. The action in a billiard scales with \^E, such that we 
take the wavenumber k := V2E/H as the integration variable and determine |p(L)| 2 from 

p(L) := / dk g(k)n osc (E(k))ex.p (ikL) exp (-tk) . (45) 
Jo 

Here the factor g(k) = H 2 k gives the measure with respect to the wavenumber. The fading function 
exp(— tk) is introduced to reduce the significance of higher eigenvalues and it has been chosen 
in order to make the analytical calculations feasable. t has to be chosen appropriately in order 
to incorporate the finiteness of the available energy range. Taking the exact eigenvalues up to 
E max = 100 000 (h = 1) calculated according to the method described in sec. [| one deals with 
n max = 35 169 levels for the ellipse parameter a = l/v2- Throughout this section we choose a = 
l/v/2 for reasons that will become clear below. The condition to be imposed on t is t 3> ln2//c max 
with A; max = y/2mE max /H. We found it adequate to set t = 0.025. 



Figure 16 shows |p(L)| 2 versus the length L. The tick marks above give the lengths of periodic 
orbit representatives of resonant tori. The spectrum shows equally spaced clusters of contributions. 
The small enclosed Figure gives a magnification of the range [4.0, 5.4] around the first cluster. 
Here the peaks correspond to type R resonant tori, with winding number l/n p , the whispering 
gallery orbits. Their lengths accumulate in the length of the sliding orbit along the billiard 
boundary Ob that can be considered as the fj, p —* oo limit of these orbits. The amplitudes of 
the contributions from these orbits decrease with growing number of reflections /i p , damped by 
(1 + /i 2 )™ 1 / 2 , the diverging curvature, and the diverging frequency, see (|33|). The other clusters 



in Fig. 16 lie around integer multiples of L^. Here multiple traversals of the whispering gallery 
orbits with winding number l/fj,„ and new ones with winding number 2//i p , 3/fi p (coprime) etc. 
accumulate. Contributions of type O tori occur only sparsely in the spectrum as compared to the 
type R tori. They have no accumulation orbit like the type R tori. The shortest orbits lying on 
a type O resonant torus have length L = 9.23 followed by L = 13.18 and L = 14.67. This means 
that the low part of the spectrum is dominated by type R tori and the clusters they produce. The 
whispering gallery orbits lead to an infinite number of (families of) periodic orbits with finite action, 
thus violating the generic growth behavior of the number of periodic orbits in integrable systems. 



In addition to resonant tori contributions Fig. 16 also shows peaks at lengths corresponding to the 



unstable orbit o u and to the stable periodic o s and multiples thereof. 
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Fi gure 16: Power spectrum |p(L)| 2 for a = i/%/2. The smaller figure shows a magnification of the range [4.0, 5.4] 
enclosed by the arrows. The tick marks above mark resonant tori, the tick marks and labels below mark isolated 
periodic orbits. 



In the following we will have a closer look at the amplitudes corresponding to isolated periodic 
orbits, and to the amplitudes of the low resonant tori. For this purpose we consider the real 
part, i.e. the cosine transform, of p(L) because it reveals much more information than the absolute 
value. Richens showed in [jl6| that as a limiting case the uniform version of the Berry- Tabor 
summation contains contributions of stable isolated periodic orbits equal to the corresponding 



terms in Gutzwiller's trace formula [34]. He also suggests that the unstable isolated orbits should 



appear in a similar way, which is rigorously shown for the ellipse billiard in [19]. Thus the qth. 
traversal of o s and o u contribute 



T(E) 



irhy/\det (Mi 



: COS 



S(E) 

h 



7T 



(46) 



to the density of states. Here T(E) is the period, S(E) the action of the periodic orbit, M 
the reduced 2x2 monodromy matrix describing its stability and [i q the Morse index of its qth 
traversal. For the calculation of Morse indices see e.g. [37]. Since for the stable orbit we have 
fig = 4q + 1 + 2[gu>], where [] denotes the integer part it is not immediately possible to factor out 
q. For o s and o u trM = 2 — 4L (k± + K2) + Aki^L 2 , where L is the flight length between two 
consecutive reflections and k\ and K2 are the curvatures at the reflection points with a positive sign 

Vl — a, 2 , while 



in case of a convex billiard boundary l38f[ . For o s we find L = L s /2 and k\ = K2 
for o u the trace is always larger than 2 because L = L u /2 = 2 and n\ = K2 = (l — a 2 ) 
winding number w G [0, 1) for o s and the stability exponent u for o u are given by 



The 



exp(±«27ra;) 
exp(±u) 
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± V(trM) 2 -4 + trM 



(47) 



where the corresponding matrix M has to be inserted. The winding number w obtained from the 
eigenvalues of the monodromy matrix M is the same as in Eq. (|9|). Now we can write the stability 
term in (|46| ) as y^det (M 9 — 1)| = 2] sm(qirw)\ or 2 sinh(qu/2), respectively. Matching the signs of 



the two sine functions for o s one can rewrite ([!(]) as 

T(E) sm{qS/h) 

2nh sm(qirw) 
T{E) cos(q(S/h- n u 7r/2)) 

2ttH sinh(gn/2) 



for o s , elliptic, 
for o u , hyperbolic, 



(48) 



where fi q = qfj, u has been used for o u , with fj, u = 4+2. It is important to notice that the amplitudes 
of the contributions of isolated periodic orbits and of resonant tori differ in the power of h, the 
former is proportional to 1/h, the latter to I//1 3 / 2 . 

A problem arises when the winding number w becomes rational, leading to a divergent amplitude 
for the contribution of the corresponding number of traversals where neighboring trajectories of o s 
are closed in phase space. In order to study this phenomenon we take a = l/\/2 in this section, 
slightly different from a = 0.7 as before. Then the winding ratio for the stable orbit Eq. fll^ ) 
becomes w = 1/2. As worked out by Richens, in this case the thin resonant torus surrounding the 
periodic orbit rather than the periodic orbit alone determines the contribution to the density of 
states. Equation fliq) then has to be replaced by 

J T(E) cos (g (S(E)/h -7TW- |qg) + fg) 

^ 2 27r(l + ™ 2 ) 3/4 ^Cp)T ' 1 ) 

The amplitude is again proportional to 1/h 3 / 2 signaling a torus contribution. 

Inserting the semiclassical results for n osc (E) into Eq. (f45|) and taking the real part we obtain 

poo 

/ dkg(k)n osc (E(k))cos(kL)exp{-tk) ~ V ^p.obj.(^) (50) 

p.obj. 

where the summation on the right hand side runs over all "periodic objects", i.e. resonant tori, 
isolated periodic orbits and thin resonant tori in cases where the isolated stable periodic orbits 
become resonant. Taking the fixed EBK phases in Eq. (33) all the different kinds of contributions 



^p.obj. can be calculated analytically. The scaling properties of the action variables in equations 
(||) and (||), the amplitudes in equations (|33j), (53), and (|4^) allow for a scaling with respect to 
the wavenumber k. Hence the semiclassical results together with ( |45l) lead to a summation over 
integrals of the form 

A(L) := A J dk cos (Zfc + n-J cos (L/c) exp (-t/c) k a (51) 

with n G Z, A, L, L, t G R, t > and cr = 0, 1/2. Defining the functions 

/ cos (L) := r(f7 + 1) ^ cos ((a + 1) arctan( )) (52) 

({L-LY + t 2 ) 2 

r(cr + i) .. . L + l 

+ CT+1 cos ((a + 1) arctan( — - — )) 



(^(L + L) 2 + ^ 

and similarly f s i n (L) by replacing cosine by sine in Eq. (|52j ) the integrals for can be solved 



analytically |23] and are listed in Tab. IV. The classical quantities Jf, u», and C are understood to 



be calculated for E = 1. Then the mapping of the parameters is given by: 

= \ 2TiqnI^^= J e {-2qaii + (5) mod 8, 1/2) (53) 

V ' \ vwlMvl ! / 

for the resonant torus contributions. The degeneracy factor e is 2 for type R tori and 1 for 
type O tori. 
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Table IV: Table of the results for A(L) 
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Figure 17: Comparison of the exact to the semiclassical result for SR(p(L)) (a = l/v2). The solid line is the cosine 
transform of the exact n osc (E), the dotted line marks the semiclassical curve. (b),(c) and (d) give magnifications of 
the length ranges [2.5, 5.4], [5.4, 10.8] and [10.8, 16.2], respectively. 



L,A,n,a) = \qL s ,-—— \ (-4qw - 2qa p + P) mod 8, 1/2 (54) 



for the thin resonant torus case where qw is an integer. 



L,A,n,a ) = (qL s ,— — -j r,-2,0 (55) 

for the stable isolated periodic orbit case with qw not an integer. 

L,A,n,a) = (qL u , 27r s [nh( qu /2) ' ^~ 2q ^ mod8 '°) ( 56 ) 

for the unstable isolated periodic orbit case with positive trace. 

The results obtained from these formulas are shown in Fig. [17]. Here resonant tori with winding 
number w = \x$j pL p and fj,^, fJL p G {1, 50} are included as far as they are realized in phase space. 
Figure (a) again shows the clustering of the contributions any time a multiple of the sliding orbit 



length Lfc is met. In Figures (b), (c) and (d) we show magnifications of the ranges between the 
(q — l)th and qth. multiple of Lb for q = 1,2,3. The small ticks are labeled by the length of the 
multiple traversals of Of,, o s and o u that can be found in that range. The large ticks labeled by 
numbers above belong to a selection of resonant tori listed in Tab. [V]. The agreement between the 
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Table V: Data for the periodic orbit representatives of resonant tori marked in Fig. [l7]. w is the winding ratio, L 
the length and q the number of traversals. 



exact and semiclassical curve is remarkably good. It gets a little worse when L becomes larger and 
the density of peaks grows. Then the sum gives an enormous mixture where the distinction of the 
individual contributions becomes more or less impossible. In Fig. 18 we show magnifications around 
some individual periodic objects. The third row shows how the amplitude of the contribution of the 
traversal of o s alternates in magnitude. For any even number of traversals the contribution is that 
of a thin torus and the amplitude is of the same order in magnitude as the torus contributions in 
the first two rows. The remaining traversals contribute as ordinary stable isolated periodic orbits. 
The fourth row shows the fast decrease of amplitudes of the unstable periodic orbit o u with a 
growing number of traversals. This feature is familiar from hyperbolic systems where all orbits 
are of this kind and the exponential decay of amplitudes with orbit length justify, e.g., the cycle 
expansion of quantum mechanical as well as classical dynamical Zeta functions (see, e.g., [p9|]). 

We performed the same calculations for a = 0.7, where o s is close to resonant and also for 
a = cos(-7r(v / 5 — l)/4) where w becomes equal to the golden mean, i.e. o s becomes as far away from 
resonant as possible. The results not represented here indicate that any time a q-fold traversal 
of o s has neighboring trajectories that are almost closed in phase space it is better to replace the 
stable orbit contribution by a thin torus contribution. For a = 0.7 this means we have a strong 
contribution for the first even traversal of the isolated stable periodic orbit which can be viewed as 
a result of the real complex orbit corresponding to w = 1/2. For a = cos(7r(v5 — l)/4) this is the 
case for q = 3,5,8,..., the denominators of the continued fraction approximations 2/3,3/5,5/8,... 
for the golden mean. 

The version of Berry- Tabor formula in Eq. (|34] ) gives a uniform expression that is valid for all 
winding numbers of stable isolated periodic orbits. In the non-rational case the contributions to 
the inverse spectrum can no longer be calculated analytically. In the inverse spectrum was 
calculated numerically. We restricted ourselves to a resonant case in order to be able to obtain 
analytical results. 



7 Conclusions and outlook 

From the classical point of view the billiard in the ellipse is a typical integrable system: the 
frequencies change from one torus to another and there exist both, stable and unstable isolated 
periodic orbits, the latter leading to a separatrix on the energy surface. Extending the classical 
mechanics to the complex plane, we introduced three kinds of complex orbits: (i) tunneling orbits 
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Figure 18: Comparison of the exact (solid) and semiclassical (dotted) result for 5ft(p(L)) in the neighborhood of the 
following periodic objects (a) (i)-(iii) traversals of the type R torus with w = 1/3, (b)(i) and (ii) traversals of type O 
torus with w = 2/3, (b) (iii) type O torus with w = 3/5, (c) traversals of the stable periodic orbit o s and (c) (i)-(iii) 
traversals of the unstable periodic orbit o u . The width for all pictures is 0.4, the heights are 140 or 20. 



and (ii) scattering orbits, both with imaginary action ("imaginary complex orbit"), and (iii) orbits 
with complex turning points, complex momentum and real action ("real complex orbits"). The 
classification of these orbits in terms of the position of the turning points and their connection by 
a Stokes or anti-Stokes line applies to all separable systems. The action of these orbits appears 
naturally in the semiclassical treatment of these systems. 

The separation of the Schrodinger equation leads to special cases of the spheroidal wave equa- 
tion corresponding to "spin" ±1/2. With a simple transformation the two coupled boundary value 
problems can be turned into a form suitable for the application of a standard shooting method 
in order to efficiently calculate the eigenvalues. The two discrete spatial symmetries of the ellipse 
lead to the distinction of four parities, while the time reversal symmetry leads to an asymptotic 
degeneracy of tori involving a rotational degree of freedom. As expected from the WKB approxima- 
tion the wave functions are concentrated on the projection of the classical torus onto configuration 
space. The isolated unstable periodic orbit induces a less pronounced "scar" for states close to the 
separatrix. Additionally these wave functions strongly localize at the focus points of the ellipse. 
An explanation of this behavior is under investigation. 

The eigenvalue spectrum can best be understood when transformed to the classical action space, 
where eigenstates are located either on simple EBK lattices far away from the separatrix or in a 
more complicated transition regime close to the separatrix. In order to semiclassically describe 
the transition between the EBK lattices in the vicinity of the separatrix a uniform quantization 
scheme was employed, which smoothes the discontinuity of the Maslov indices at the separatrix. 
The imaginary complex orbits are used in this context to define a smooth Maslov phase, replacing 
the discontinuous Maslov index. The semiclassical states are obtained as the intersections of a set of 
lines in classical action space. Numerically these states agree quite well with the exact results, even 
for low quantum numbers. The relative error shows regular patterns, except for separatrix states, 
and it tends to be large for states with at least one small quantum number and for states close 
to the separatrix. A state with given quantum number is located inside a "semiclassical quantum 
cell", and specifying its parity narrows its position down to a "parity box" of width h/2. This 
interpretation of the uniform quantization allows for an illuminating picture of the situation for 
states close to the separatrix. Without uniformization, parity boxes intersected by the separatrix 
may contain a wrong number of EBK states. 

The Berry- Tabor trace formula rests on EBK quantization. In the original uniform version it 
contains contributions from real complex orbits, but imaginary complex orbits are not incorpo- 
rated. A treatment ignoring the separatrix produces faulty peaks in the semiclassical density of 
states. They occur for states that belong to parity boxes with the wrong number of EBK states. 
Incorporating the imaginary complex orbit by a smoothly varying Maslov phase leads to a modified 
trace formula. In the semiuniform WKBBT version these corrections are ignored in the stationary 
phase approximation, resulting in a sum over resonant tori as in the original version. This, however, 
cannot correct all the peaks from states belonging to parity boxes with the wrong number of EBK 
states. In order to obtain satisfactory peaks for these states we had to fully take into account the 
non-constant Maslov phase. This leads to a uniform WKBBT trace formula quite similar to the 
original one, with a changed effective energy surface now depending on the energy. The effective 
resonant tori of this effective energy surface do not correspond to periodic orbits of the billiard. 
With this sum over (in general) non-periodic orbits we were able to produce the correct peaks for 
all states. The numerical effort increases considerably because the stationary points have to be 
found anew for every energy. Even though the uniform WKBBT trace formula might not be a 
useful tool for the practical calculation of eigenvalues, it can give some hints on how to incorporate 
complex orbits into formulas of this type. It remains an open question how a uniformization should 
come about in terms of a Green's function derivation of the Berry- Tabor trace formula jSj]. It is 
by no means obvious how the variable Maslov phase should be emulated by, say, a summation over 
complex orbits. 

In the study of the length spectrum of our integrable system we focused on the case where 
the stable orbit has a rational winding number. In this case the divergent Gutzwiller term has 



to be replaced by a thin torus contribution obtained by Richens. The four types of dominant 
contributions in the trace formula can be analytically Fourier transformed, such that the inverse 
spectrum can be explicitly written as a sum over periodic objects. We found a remarkably good 
agreement, even though in this approach only the real complex tori are partially incorporated in 
the Richens term. It is an open question why the inverse spectrum is so much less sensitive to the 
presence of the separatrix in our case. 

Extensions of our work can be done in two directions. On the one hand one can take the elliptic 
billiard as a starting point for the penetration into the nonintegrable regime by a deformation of 
the billiard boundary. This is the content of [IS]. On the other hand one can consider the three 
dimensional billiards in the ellipsoid [10, 41 1, starting with the cases of prolate and oblate spheroids. 
A first step in this direction can be found in p^]. The Berry- Tabor trace formula for a system with 
three degrees of freedom is much more involved, because the resonant tori can no longer be labeled 
by a single rational winding number. Instead one has to find a complicated set of resonances 
on the energy surface [41]. For these 3D billiards, the exact quantum mechanical spectrum, its 
semiclassical quantum cells, and and its length spectrum are currently under investigation. 
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